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Abstract 

A few per cent fraction of antiprotons stopped in helium survives for an enormous 
time (up to tens of microseconds) in comparison with the usual lifetime (10~^'^s) 
of these particles in matter. The explanation is that antiprotons are captured in 
the metastable antiprotonic helium atoms ^''^Hepe carrying an extremely large total 
angular momentum L ~ 30 — 40. Initial populations, level lifetimes and very precise 
values of the transition energies were obtained in experiments on the resonant laser- 
induced annihilation. Analogous long-lived systems were observed in experiments 
with negative kaons and pions. 

The purpose of this report is to present the results of calculation of eigenenergies, 
radiative transition rates, energy-level splitting due to relativistic interactions and 
Auger decay rates of antiprotonic helium atoms within the variational approach. 



1 Introduction 



The lifetime of negatively charged heavy particles (such as tt~,K~,P) stopped in matter 
would be about 10~^^s due to nuclear absorption after exotic-atom formation. Recently 
it has been found that a small fraction (a few per cent) of kaons pions and an- 
tiprotons stopped in helium media survive for an enormous time (up to tens of 

microseconds for antiprotons). The reason for this unusual behavior is the formation of 
a metastable three-body system consisting a helium nucleus, an electron and a negative 
heavy particle. The prediction of these systems were made as early as thirty years ago 
by Condo ^ to explain a considerable difference of pion lifetimes in liquid helium and 
hydrogen. Later on, this qualitative prediction was supported by the variational calcula- 
tions of Russell 10. One can mention also the first calculation of the antiprotonic helium 
atom using the Born-Oppenheimer approximation [0. Since antiprotons are stable and 
lifetime of an antiprotonic helium is the largest among these exotic systems, the exper- 
imental data concern mainly the antiprotonic helium atom. The modern experiments 
on the resonant laser-induced annihilation P, ^ have initiated thorough investigations 
of these unusual systems. Initial populations, level lifetimes, and very precise values of 
the transition energies (with a relative accuracy < 10~^) have been obtained using this 
remarkable method. 

The antiprotonic helium can be considered as a counterpart of the usual helium atom 
with one electron replaced by an antiproton. A large angular momentum L ~ (fx/m^y^'^ 
{li is the reduced mass of the antiproton-nucleus subsystem and ms is the electron mass) 
provides that the antiproton-nucleus and electron-nucleus distances are approximately 
equal. At the same time, this system can be considered as an exotic diatomic molecule, 
where one nucleus is negatively charged [|T0|] . 

While there is a number of processes destroying the antiprotonic helium atom, the 
following considerations explain its very long lifetime. The most part of an extremely 
large total angular momentum L ~ 35 — 40 belongs to the pair of heavy particles and 
annihilation of an antiproton is inhibited by a large centrifugal barrier. The Auger decay 
is inhibited due to sufficiently large values (Aq > 4) of the angular momentum of an 
outgoing electron. Usual mechanism of the de-excitation by the Stark mixing is not 
appropriate for the three-body system due to the lack of degeneracy. The coUisional de- 
excitation by surrounding He atoms is suppressed due to the screening of an antiproton 
by the electron in an antiprotonic helium. The only remaining de-excitation mechanism 
is multistep dipole radiative transitions, whose rates of each step are of order yus~^. The 
discussion of theoretical calculations on antiprotonic helium atoms and related topics can 



be found in 11 



After the formation of the antiprotonic helium its time evolution is determined by 
the processes of radiative transition and Auger decay. The precise description of energy 
spectra requires to take into account the minor effect of relativistic interactions. The 
spin-dependent part of relativistic interactions gives rise to splitting of energy levels and 
sufficiently large splitting values can be measured experimentally. The direct variational 
method provides a possibility in the framework of one approach to calculate eigenenergies, 
energy-level splitting due to relativistic interactions, radiative transition rates and Auger 
decay rates of antiprotonic helium atoms. The results of these calculations are presented 
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and discussed in this report. The isotopic effect is considered by comparison of the 
calculated properties of ^Hepe and ^Hepe. 



2 Variational method 

The nonrelativistic Hamiltonian of the antiprotonic helium atom in the triangular coor- 
dinates is 

11 1 2 2 1 

H = - — Ar - — Ap Vr ■ Vp + I (1) 

Zjj, Ijii nil r p |r — p| 

where l//i = 1/mi + l/m2, l//ii = 1/mi + l/ms and mi, ri, m2, r2, ma, are masses 
and coordinates of the helium nucleus, antiproton and electron, respectively. Apart from 
the triangular coordinates r = r2 — ri,p = — ri, the Jacobi coordinates r = r2 — 
ri, Pi = rs — (miFi — m2r2)/(mi + m2), corresponding momenta p = — iVr, q = — zVp, 
qi = —iVp^ and angular momenta 1 = [rp], A = [pq], Ai = [Piqi] will be used to 
simplify the notation. 

Since antiprotonic helium atoms are unstable against the decay to the Hep + e channel, 
the variational method cannot be directly applied to the calculation of energies and wave 



functions. Therefore, according to |12| the approximate Hamiltonian H^n = Pl^HPln 
which is the projection of the hamiltonian H onto the closed-channel subspace was con- 
structed to calculate the N-th energy level Eln of the system for the total angular momen- 
tum L. Explicitly Pln is constructed as a projector onto the subspace of eigenf unctions 
of relative angular momenta 1, A (or 1, Ai), belonging to a limited set of /, A (or /, Ai) 
eigenvalues. While using the Jacobi coordinates this set is defined by the condition I > Iq, 
where Iq is the largest pair angular momentum satisfying the inequality ei^ < Elm and 

is the energy of the hydrogen-like ion ^'^Hep with the angular momentum /q- If the 
triangular coordinates are used, this set can be defined by the condition A < Aq, where Aq 
is the smallest angular momentum satisfying the inequality ^l-Ao < -^LAf- By definition, 
Aq is also a multipolarity of the Auger transition, i. e. the smallest angular momentum 
of the outgoing electron. These conditions describe the natural way to provide the ap- 
proximate Hamiltonian Hl^ to have at least N eigenvalues below the boundary of the 
continuous spectrum and the variational method can be applied to solve the eigenvalue 
problem for lowest states of Hi^. 

It is worthwhile to mention that the / = L projection is the main part of the wave 
function and the contribution oil ^ L projection rapidly decreases with increasing |/ — 
For this reason the approximate Hamiltonian Hl^ provides an accurate calculation of 
energies El^ and wave functions ln for an metastable states (Aq > 4) of antiprotonic 
helium atoms. In fact, the described conditions on /, A are more restrictive than necessary 
for an application of the variational method. Practically the accuracy of calculation can 
be improved taking into account A > Aq (or I <Iq m the case of the Jacobi coordinates) 
components of the wave function. While using these components a convergence of the 
calculated values with increasing the number of trial functions will be investigated to 
provide the reliable results. 
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3 Energy levels 



The variational method described in the previous section was apphed to determine eigen- 
functions and eigenenergies of the equation 



(Hln — Eln)'^ LN = 0. 

A set of simple variational trial functions of the form 

3^,i^^(f , py+'p^expi-a^r - hp), 



XnklXii^i P) 



(3) 



where yj^^ir^p) are bispherical harmonics of angular variables, was used in the calcu- 
lations. It is essential that yl^^\r,p) are eigenfunctions of 1, A and this form of trial 
functions allows the easy application of the projection method described in Section 0. 

Up to 1000 trial functions containing up to 25 bispherical harmonics were used 
in the calculations. A set of nonlinear parameters a„, for simplicity was chosen in 
the form a„ = aoa"',bk = h^P^ and the variation of parameters ao,a,6o,/? was used to 
minimize energy values. As it is usual for the variational method, the precision falls down 
with increasing the excitation number N . To reach more precise energy values, a set of 
trial functions was limited by the condition / > L — Aq or even less restrictive condition 
/ > L — Ao — 1 instead of A < Aq. Convergence of the calculated energies with increasing 
the number of trial functions was obtained and results are presented in Table |I]. Some 



Table 1: Calculated energies ELj^ioj.u. 
in the range 31 < L < 37 



for five lowest energy levels of the ^'^Hepe systems 



'^Hepe 


L 


Eli 


El2 


El3 


El4: 


El5 


32 


-3.3534091 


-3.2272500 


-3.1161769 


-3.0183048 




33 


-3.2158963 


-3.1049859 


-3.0075402 


-2.9219136 


-2.8355338 


34 


-3.0931201 


-2.9959130 


-2.9107010 


-2.8359650 


-2.7685835 


35 


-2.9836255 


-2.8988545 


-2.8247010 


-2.7597140 


-2.7007897 


36 


-2.8862911 


-2.8126761 


-2.7483976 


-2.6920660 


-2.6389254 


37 


-2.7999636 


-2.7364076 


-2.6809542 


-2.6320110 


-2.5853018 


'^Hepe 


L 


Eli 


El2 


El3 


El4 


El5 


31 


-3.3484555 


-3.2190748 


-3.1056452 


-3.0061026 




32 


-3.2073388 


-3.0940418 


-2.9949034 


-2.9082665 


-2.8309327 


33 


-3.0817523 


-2.9829571 


-2.8967378 


-2.8214095 


-2.7539447 


34 


-2.9702506 


-2.8844813 


-2.8097982 


-2.7445651 


-2.7007897 


35 


-2.8714926 


-2.7974268 


-2.7330406 


-2.6763097 


-2.6389254 


36 


-2.7843183 


-2.7207232 


-2.6654633 


-2.6161392 


-2.5853018 



information on this calculation can be found also in |12|. Energies of the L,N states of 
^Hepe and L — 1,N states of ^Hepe are very close to each other. 
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By using the method of laser- induced resonant annihilation, till now, wavelengths of 
five transitions (35,4 ^ 34,4) §, (34,3 ^ 33,3) (35,2 34,4), (34,3 ^ 33,5), 
(35,3 ^ 34,3) [|1|] in ^Hepe and two transitions (34,4 ^ 33,4), (33,3 ^ 32,3) H in 
^Hepe were measured with high accuracy. The calculated and experimental wavelengths 
are in agreement with an accuracy not worse than 5 ■ 10~^ and more elaborate trial func- 
tions are needed to reach a higher accuracy in the energy calculation. Recently, V. I. Ko- 
robov has obtained more precise spectra and transition wavelengths of an antiprotonic 
helium, using the correlated trial functions in variational calculations []T5|. 



4 Radiative transitions 

Due to large lifetimes of metastable states against the Auger decay and collisional de- 
excitation the radiative transitions become the most important to describe the evalution 
of the antiprotonic helium atom. As far as only dipole transitions are significant, the 
total angular momentum changes by unity in each transition. Thus, the system looses 
the angular momentum and energy step-by-step and finally reaches the state with a large 
Auger decay rate. 

The rate of the dipole transition LN LiNi is given by 

w = 3(2L^i) l^(-^^^ - EL,NAf\Md\^^s-\ (4) 

where a = /he is the fine structure constant and the reduced matrix element of the 
dipole transition is defined in the form 

Md= {^L^NA\r + pW^ln). (5) 

Variational wave functions ln and energies Eln have been obtained as described in the 
previous section and used in (|), (|]) to calculate the radiative transition rates. These 
results are presented in Table ^. A convergence of the calculated values with increasing 
the number of trial functions provides an estimate of the relative accuracy on the level of 
one per cent. The radiative transition rates calculated in papers |T^, [jl6|, [l^ are fairly 
close to each other and to present results. The important feature is the predominance 
of transitions between states of the same A^, i. e. the approximate conservation of the 
excitation number in the radiative transitions. Thus, radiative cascades in an antiprotonic 
helium proceed almost independently along the chains of states of fixed A^. However, 
for higher excitation numbers the probabilities of interchain transitions become more 
significant and such transitions will be taken into account for the cascade description. 

As in the energy calculation, there is a correspondence of the L, A^ state of ^Hepe 
and the L — 1,N state of ^Hepe and the radiative transition rates of the corresponding 
states are almost the same. At the same time, the ^Hepe transition rates systematically 
exceed those of ^Hepe in accordance with the difference of experimental lifetimes of the 
antiprotonic helium atoms [Q. However, the explanation of this difference is not simple 
because the mean lifetime depends also on another processes and a significant difference in 
populations of the corresponding states was found in the recent experiments for different 
isotopes fl^ . 
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Table 2: Radiative transition rates w(10^s ^) from the LN to LiNi state of the ^'^Hepe 

systems. Only transition rates w > 10^s~^ are presented. 

^Hepe ^Hepe 

transition w transition w transition w transition w 

37,1^36,1 6.733 35,5^34,5 3.143 36,1^35,1 7.255 34,5^33,4 1.036 

37,2^36,1 0.159 34,1^33,1 7.744 36,2^35,1 0.187 34,5^33,5 4.532 

37,2^36,2 5.859 34,2^33,2 7.554 36,2^35,2 6.255 33,1^32,1 8.483 

37,3^36,2 0.441 34,3^33,1 0.206 36,3^35,2 0.511 33,2^32,2 8.240 

37,3^36,3 4.907 34,3^33,3 6.926 36,3^35,3 5.182 33,3^32,1 0.301 

37,4^36,3 0.742 34,4^33,1 0.288 36,4^35,3 0.887 33,3^32,3 7.539 

37,4^36,4 3.990 34,4^33,2 0.137 36,4^35,4 4.135 33,4^32,1 0.476 

37,5^36,4 1.218 34,4^33,3 0.355 36,5^35,4 1.642 33,4^32,2 0.205 

37.5^36.5 2.917 3 L 1^33.1 6.073 36.5^35.5 2.931 33 32.3 0.387 

36,1^35,1 7.249 34,5^33,1 0.355 35,1^34,1 7.846 33,4^32,4 6.577 

36,2^35,2 6.531 34,5-^33,2 0.365 35,2^34,1 0.102 33,5^32,1 0.748 

36,3^35,2 0.333 34,5^33,4 0.728 35,2^34,2 7.019 33,5^32,2 0.571 

36,3^35,3 5.644 34,5^33,5 4.840 35,3^34,2 0.384 33,5^32,4 0.807 

36,4^35,3 0.637 33,1^32,1 7.720 35,3^34,3 5.995 

36,4^35,4 4.685 33,2^32,2 7.859 35,4^34,3 0.747 

36,5-^35,4 1.077 33,3^32,1 0.648 35,4^34,4 4.937 

36,5->35,5 3.323 33,3^32,3 7.492 35,5->34,4 1.304 

35,1^34,1 7.629 33,4^32,1 1.166 35,5^34,5 3.666 

35,2^34,2 7.123 33,4^32,2 0.570 34,1^33,1 8.276 

35,3^34,2 0.205 33,4^32,3 0.151 34,2^33,2 7.700 

35,3^34,3 6.331 33,4^32,4 6.798 34,3^33,2 0.236 

35,4^34,1 0.137 33,5^32,1 3.503 34,3^33,3 6.797 

35,4^34,3 0.509 33,5^32,2 1.537 34,4^33,1 0.169 

35,4^34,4 5.327 33,5^32,3 0.321 34,4^33,3 0.582 

35,5^34,1 0.168 33,5-^32,4 0.434 34,4^33,4 5.762 

35,5^34,2 0.187 33,5^32,5 0.730 34,5^33,1 0.199 

35,5^34,4 1.059 34,5^33,2 0.244 
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Lifetimes of the (L, A^) = (35, 4) and (L, N) = (34, 3) states of ^Hepe system were 
determined experimentally by the method of resonant laser-induced annihilation , . 
Table ^ contains the experimental decay rates and theoretical radiative transition rates 
of these states. One of the most probable reasons for the difference of experimental and 



Table 3: Experimental decay rate and theoretical radiative transition rates (lO^s ^) for 
two states of the '^Hepe system 



L,N 


Experiment 


m 


m 


present 


35,4 


0.72 ±0.02 @ 


0.614 


0.619 


0.597 


34,3 


1.18 ±0.04 ll 


0.734 


0.754 


0.713 



all the theoretical values indicates the substantial contribution of additional nonradiative 
decay channels. As a density dependence of the lifetime of the (34,3) state in ^Hepe has 
been found in recent experiments [|14|, collisions with surrounding atoms contribute to 
the decay processes. An extrapolation of the new experimental data to zero density gives 
0.91 ± 0.15yus~^, which is close to the theoretical radiative transition rate. The density 
effect probably takes place also for the 35,4 state and the remaining discrepancy may be 
caused by an additional decay process, e. g. the Auger decay. 



5 Energy— level splitting due to relativistic interac- 
tions 

The precise measurement of transition energies of antiprotonic helium atoms in recent 
experiments on the laser- induced resonant annihilation [Q, invokes the theoretical 
description of energy spectra with a comparable accuracy. That description of energy 
spectra requires to take into account the relativistic corrections of an order of to the 
pure Coulomb interaction. 

Since the contribution to energies from relativistic interactions depends on the antipro- 
ton mass, charge and magnetic moment, these calculations can be used for the precise 
determination of the antiproton properties. This knowledge is essential in testing the 
fundamental symmetry principles |18]. 



The spin-dependent part of the relativistic interactions gives rise to splitting of energy 
levels, and each single transition turns into a multiplet. Sufficiently large distances be- 
tween lines in the multiplet can be measured experimentally. It is worthwhile to mention 
that the resolution in current experiments is about lOGHz and without much difficulty 
can be improved to IGHz [^]. As it will be discussed below, due to the interaction with 
electron spin, antiprotonic helium energy levels split into two multiplets and the interac- 
tion with nuclei spins provides a minor splitting within each multiplet. Calculation of the 
former large splitting is discussed in this section. More details of this calculation were 
presented in . 

For each pair of particles i,j in the three-body system the relativistic correction of 
an order of to the pure Coulomb two-body potential can be described by the Breit 
interaction U^P. The correction to the kinetic energy of an order of for each particle 
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I IS 

AT, = -^4 (6) 

8 mf 

Full relativistic correction of an order of to the three-body nonrelativistic Hamil- 
tonian is a sum of Ulf^ for all pairs of particles and ATj for all particles 



i>3 



In the definition of the expressions U^f^ , ATj in eq. (|^) it is proposed that particles 
momenta pj will be taken in the center of mass system of three particles . 

The interaction given in (|^) conserves the sum J = L + Sj of the total angular 
momentum L = 1 + A and particle spins Sj. Each level of the nonrelativistic Hamiltonian 
splits into four and eight sublevels for ^Hepe and ^Hepe systems, respectively. Due to very 
small mass ratios 1713/1711,7713/7712, the largest contribution to the energy splitting comes 
from the interaction with the electron spin S3. Taking into consideration only terms 
responsible for the splitting in (|^), this part of relativistic interaction can be written as 
follows: 

Hs = a\\xs3+ ^ [r-p,q]s3- (8) 

1 2 

[r - P, p]s3 + ^[p, p]S3) 



m2|r — p\-^ 77iip-^ 

While the last two terms in (j^) are inversely proportional to the masses of heavy par- 
ticle mi^2, their contribution to the energy splitting is nevertheless comparable to the 
contribution from the first two terms for the following reasons. The small mass factor is 
compensated in part due to the large angular momentum / ~ L of heavy particles. At 
the same time, only small components of the wave function corresponding to the nonzero 
electron angular momenta A 7^ lead to a nonzero splitting value from the first two terms 
in ®. 

The interaction Hg, given in (|^), conserves the sum j = L + S3 of the total angular 
momentum L and electron spin S3 and splits each level into two sublevels, corresponding 
to the eigenvalues j = L ±1/2. The part of the interaction depending on heavy particle 
spins removes the remaining degeneracy and splits each j = L ± 1/2 sub level further into 
two or four levels for the '^Hepe and ^Hepe systems, respectively. Values of this secondary 
splitting are much smaller in comparison with the splitting arisen due to the interaction 
with the electron spin 

Since the splitting is small in comparison with energy differences between states of 
different L values, the energy shift A^t^at for the state with quantum numbers j, L, N can 
be found in the first order of perturbation theory in Hg 

AjLN= {^jLnIHsI'^jln), (9) 

where the wave function "^jiN of the j, L, N state is the vector production of the nonrel- 
ativistic wave function "^ln and spin function describing the dependence of the electron 
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spin. Level splitting AE^n = AL+i/2LAf ~ ^l-i/2LN is a difference of shifts (|) for 
j=L± 1/2. 

Due to smallness of the relativistic interaction, radiative transitions proceed only be- 
tween states of the same j. For this reason, each spectral line of the transition from 
the state LiNi to state LjNj is to be split into a doublet with the interline distance 

Ap = AEl^n,- AEL^Nr 

Splitting values AE^q for a number of states of the ^''^Hepe systems in the range of 
experimentally observed values of the total angular momentum L have been calculated 
as described above by using variational nonrelativistic wave functions "^ln- These values 
are presented in Table 

Table 4: Splitting values AEln (10~^au) of the lowest levels in the ^''^Hepe systems. 



'^Hepe 


N 


L=32 


L=33 


L=34 


L=35 


L=36 


L=37 


1 


-1.10 


-1.15 


-1.15 


-1.14 


-1.12 


-1.09 


2 


-1.12 


-1.09 


-1.08 


-1.07 


-1.04 


-1.00 


3 


-1.01 


-1.02 


-1.00 


-0.98 


-0.94 


-0.90 


4 




-0.94 


-0.94 


-0.90 


-0.86 


-0.82 


5 






-0.93 


-0.90 


-0.84 


-0.81 


^Hepe 


N 


L=31 


L=32 


L=33 


L=34 


L=35 


L=36 


1 


-1.20 


-1.16 


-1.19 


-1.19 


-1.18 


-1.14 


2 


-1.14 


-1.19 


-1.15 


-1.12 


-1.08 


-1.04 


3 


-1.08 


-1.06 


-1.05 


-1.04 


-1.00 


-0.98 


4 




-0.97 


-0.92 


-0.86 


-0.85 


-0.81 



As it follows from expression (^, the form of the wave function at small interparti- 
cle distances is the most important in evaluating the integral (^. Convergence of the 
calculated splitting values AE^j^ provides a few per cent relative accuracy. It is worth- 
while to mention that due to a better description of the small states in the variational 
method, the accuracy of calculation for the L, N state increases with decreasing A^. At 
the same time an accuracy of calculation decreases with decreasing L due to decreasing 
the multipolarity Aq of the Auger decay for large L. Moreover, it is impossible to trace 
the convergence in the case of short-lived states due to a small multipolarity Aq < 3 of 
the Auger decay. This problem is closely connected with a large natural width of these 
states which exceeds significantly a splitting value. Also, the variational procedure meets 
some difficulties in describing the short range behavior of the wave function for states 
with large enough A^, especially, in the ^Hepe system. These are the reasons to omit the 
above-mentioned cases in Table |^. 

The last two terms in eq. (|^) describe the interaction of the electron magnetic moment 
with the magnetic field of heavy particles. These terms give rise to the largest contribution 
to the energy-level splitting. For a better understanding of the splitting dependence on 
L, N, this contribution is presented in Table ^ for the ^Hepe system. The contribution to 
the energy-level splitting from the first two terms in (H) is of opposite sign and smaller in 
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magnitude. Nevertheless, this contribution decreases with increasing L and compensates 
the L dependence of the last two terms in eq. (|^) providing a very slow dependence of the 
total sphtting AE^jy on L. Due to almost exact conservation of the j value in the radiative 



Table 5: Contribution of the last two terms in A^^ to the energy-level splitting AE^j^ 
(10~^au) in the ^Hepe system. 



N 


L=32 


L=33 


L=34 


L=35 


L=36 


L=37 


1 


-1.41 


-1.43 


-1.40 


-1.37 


-1.34 


-1.28 


2 


-1.39 


-1.34 


-1.30 


-1.27 


-1.22 


-1.16 


3 


-1.26 


-1.24 


-1.20 


-1.15 


-1.10 


-1.04 


4 




-1.14 


-1.11 


-1.06 


-1.00 


-0.94 


5 






-1.10 


-1.06 


-0.98 


-0.95 



transition the spectral line splitting will be found as a difference of AE^sf presented in 
Table |[ Most appropriate for the experimental measurement are the favoured transitions 
between states of the same A^, which have the largest radiative rates |]12[, |T^, [IIO 



However, the calculated splitting values are almost independent of L for a given A^, and it 
is not plausible to resolve such a small difference in splitting for the favoured transitions. 
For this reason, the experimental proposal for the nearest future |jl9[ is aimed at searching 
for the splitting in unfavoured transitions (L, N) {L — 1,N + 2). 

To measure splitting in experiments on the laser-induced resonant annihilation, the 
initial state should be long-lived. This is provided by the condition that the multipolarity 
of the Auger decay for this state is Aq = 4. The next condition is that the natural width of 
the short-lived final state will be smaller than the splitting value, and the multipolarity 
of the Auger decay for this state will be Aq = 3. The spectral line splitting Au for a 
number of suitable transitions is presented in Table |. These values are of an order of the 



Table 6: Spectral line splitting Au 
LfNf in the ^''^Hepe systems. 



AEi-j^f. —AEl n (GHz) for the transitions LiNi 



'^Hepe 


'^Hepe 


LiNi^LfNf Au 


LiNi^LfNf Au 


33,1-^32,3 -0.92 


32,1^31,3 -0.53 


34,1^33,3 -0.86 


33,1^32,3 -0.86 


34,2^33,4 -0.91 


33,2^32,4 -1.22 


35,2^34,4 -0.87 


34,2^33,4 -1.35 


35,3^34,5 -0.34 





experimentally measurable value ~lGHz. One of the recently discovered |T^ unfavoured 
transitions (35, 2 — > 34, 4) in ^Hepe is a good candidate for the splitting measurement. 

Energy-level splitting AEln decreases with increasing the excitation number A^ and an 
obvious reason for this dependence is decreasing of the wave function at small interparticle 
distances for excited states. One can mention an appreciable difference in the AElm 
dependence on A^ for the ^Hepe and ^Hepe systems. As a consequence, an appreciable 
isotopic effect appears also for the spectral line splitting Au. 
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The following considerations can be used to understand qualitatively the L, N- 
dependence of the energy-level splitting. The contribution to splitting from the inter- 
action of the electron magnetic moment with the magnetic field of heavy particles is 
described by the last two terms in the splitting interaction Hg (|^). This contribution is 
proportional to the relative momentum of heavy particles p. One can consider that the 
motion of heavy particles is approximately the same as in a hydrogen-like atom and the 
momentum p is inversely proportional to the angular momentum L. This is the reason for 
increasing this contribution with decreasing L, as presented in Table ^. The contribution 
from the first two terms in the splitting interaction Hs is connected with the electron 
rotation and is proportional to a small component of the wave function arising due to 
polarization of an electron by an antiproton. With decreasing L the antiproton moves to 
the region of increasing electron density and the polarization increases. In this way, con- 
tributions to the energy-level splitting from the last two terms in Hg and the remaining 
part of splitting interaction are of opposite sign and level off the dependence of the total 
splitting AEln on L. 

One can consider quasi-classically that the antiproton orbit became more stretched 
with increasing at a fixed total angular momentum. For this reason all the terms of 
the splitting interaction Hg decrease with increasing L and provide the N dependence 
presented in Tables ||. 



6 Auger decay 

In addition to the radiative transitions, the important decay channel of antiprotonic he- 
lium atoms is the emission of an electron with the ^'^Hep hydrogen-like ion formation, 
i. e. the Auger decay. The main feature of the Auger decay rates of antiprotonic he- 
lium atoms is the essential dependence on the transition multipolarity, i. e. the smallest 
angular momentum of the outgoing electron Aq. As discussed in Section ||, calculated 



eigenenergies E^n 10 unambiguously determine the transition multipolarities 

Ao by the condition El-Xq < Elm- 

Apart from the early calculation of Russell [0] the only progress in the calculation of 



the Auger decay rates is due to the paper [|16|. It was found that the Auger decay rate 



decreases about three orders of magnitude with increasing the smallest angular momentum 
of the outgoing electron Aq by one. 

As the Auger lifetime of states with Aq = 3 and Aq = 4 is of an order of lO^^s and lO^^s, 



respectively, and the radiative lifetime is about 10" s |I0|, [|T7|, [0, the resonant laser- 
induced transitions between states of Aq = 3 and 4 became observable in the experiments 
1^, 1^. For this reason, the most important thing is to calculate the Auger decay rates 
for the states of ^'^Hepe systems with multipolarities Aq = 3, 4. The method and results 
of the calculation are presented in and will be described in this section. One should 
mention that there are no other calculations of the Auger decay rate for the ^Hepe system 
and the presented results allow one to consider the isotopic effect in this process. 

The Auger decay rates of the antiprotonic helium atom are very small in the atomic 
scale and the very fine details of the discrete state wave functions "^^ln and continuum wave 
functions will be accurately determined. In particular, since "^ln and \l/c belong to the 
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closed- and open-channel subspaces, respectively, the orthogonality condition (\l/iAr|\l/c) = 
will be strictly fulfilled. As described in Section 0, "^ln is an eigenfunction of the 
approximate Hamiltonian H^jy and projector Pln- As a consequence, \l/c should be an 
eigenfunction of the projector 1 — Pln and can be found as a solution of the equation 

(l-PLJv)^(l-PLiv)^c = ^LJV^c. (10) 

Thus, the orthogonality requirement is fulfilled and the perturbation theory can be applied 
due to smallness of the decay rates. 

According to the Feshbach orthogonal projection method the decay rate A is 

A = , ^ |M,.r^s-\ (11) 

where l/yUa = l/(mi -|- 1712) + and the transition matrix element is 

Mtr = (^c|(l - Pln)HPlm\^ln). (12) 

Since the continuum wave function is naturally described in the Jacobi coordinates r, p^, 
these coordinates will be used in this calculation. 

The contribution to Mtr is negligible for the \l/c components, corresponding to (/, Ai) 
eigenvalues of angular momenta 1, Ai if / 7^ Iq, Ai 7^ Aq. For this reason, only the (/q, Aq) 
component will be taken into account in the calculation of \l/c. 

Due to the large centrifugal barrier the continuum wave function \E'c can be taken as a 
product of the normalized wave function of ^'"^Hep and the function /(pi) describing the 
relative motion of an electron and ^''^Hep. As the total angular momentum of the system 
is L and its projection is M, \l/c takes the form 

*,(r, pi) = A ytXa, Pi) r'o e-^^ f{p,) = <l>(r, Pi)/(pi). (13) 

In the framework of this approach the interaction of the outgoing electron and the re- 
maining ^'"^Hep ion is described by the folding potential 

Voipi) = ($1(1 - Pln)H{1 - Pl7v)|$), (14) 

and the radial function /(pi) obeys the ordinary differential equation. 

Due to the large centrifugal barrier for the outgoing electron, the Auger decay rate 
is mainly determined by the wave functions in the large pi range, where the interaction 
of the electron and ^'^Hep system is nearly a sum of the Coulomb and centrifugal po- 
tential. Indeed, the replacement of the folding potential by the Coulomb one gives rise 
to a minor change in the matrix element. This fact supports the applicability of the 



approximation (|1^). 

500-650 trial functions Xntixii^^Pi) ® variables r, p^^ have been used in the 

variational calculation of '^ln.Eln- As described above, to fulfil the strict requirement 
of orthogonality, only trial functions satisfying the condition I > Iq have been taken 
into account. The largest contributions to the transition matrix element M^^ come from 
the {l,L — I) components of the wave function and they have been treated with special 
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care. An additional set of 70-130 trial functions with specific nonlinear parameters have 
been used to describe these components. As the calculated energy values are in good 
agreement with the results of , one can conclude that the present calculations provide 
the accurate description of antiprotonic hehum atoms. Using "^ln, Eln and \l/c(r) in the 
form (0) the Auger decay rates A have been calculated according to equations (|lT]) , (|12D . 
These results are presented and compared with 16 1 in Table ^ To understand the role 



Table 7: Multipolarities Aq and Auger decay rates A (s ^) of ^'^Hepe 



^Hepe ^Hepe 



L,N Ao A A g 


L,N Ao A 


34,4 3 7-10^ 8.5-10^ 


33,4 3 2-108 


33,3 3 1.2 ■ 10« 1.5 ■ 10«[] 


32,3 3 3.4 - 10" 


32,2 3 1.1 ■ 10** 4.0 ■ 10^* 


31,2 3 3.3 -10^ 


35,4 4 - 5-10^ 6.7-102 


34,4 4 ~ 10^ 


34,3 4 ~ 1 ■ 10^ 2.7 ■ 10^ 


33,3 4 ~ 10^ 



*K. Ohtsuki, private communication 



of the wave function structure in this calculation, the largest contributions of the (/, Ai) 
components of "^ln to the transition matrix element are presented in Table ^. It is 
worthwhile to mention that the exponentially small part of the wave function in the large 
Pi region is important in evaluating the integral (]T^. As it follows from the numerical 



Table 8: Normalized contributions to the transition matrix element Mtr{l,Xi) from the 
(/, Ai) components of the wave function for two states of the ^Hepe system. 



(L,N). 


= (34,3) Ao=4 


(L,N). 


= (34,4) Ao=3 




Mt,(/,Ai) 


(/,Ai) 


Mtr{l,\l) 


34,0 


-0.28 


34,0 


0.24 


33,1 


0.51 


33,1 


-0.46 


32,2 


-0.38 


32,2 


1.22 


31,3 


1.14 





results, the essential point in the calculation of the Auger decay rates is to determine very 
fine features of the wave function. Really, the largest contribution to the transition matrix 
element comes from the smallest component corresponding to the largest possible Ai value 
(Table ||) and contributions to the Mtr from the (/, Ai) components of "ifLN compensate 
each other due to alternation of signs. Moreover, the most important in the calculation 
of Mtr is the large pi region of the configuration space, where these components decrease 
exponentially. 

As a result, uncertainty in the calculated decay rates for the multipolarity Aq = 3 
is of an order of 10 per cent. As is indicated in Table |^, the decay rates for the ^Hepe 
system in this case are in agreement with the results of [T^ . The calculated decay rate of 
the (34,4) state of this system is in fairly good agreement with experimental lifetime r= 
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15ns [§. At the same time, for the (33,3) state the calculated decay rate twice exceeds 
the experimental value (r=16.6ns) |T3[. 

On the contrary, the case of the transition multipolarity Aq = 4 is rather complicated 
for calculation due to much smaller values of the transition matrix element, and only an 
estimate of the decay rate can be obtained. These values exceed significantly the results 



obtained in 16 



The calculated decay rates of the ^Hepe system reveal the substantial isotopic depen- 
dence. In fact, decay rates of the (33,4) and (32,3) states of ^Hepe are about three times 
as large as those of analogous (34,4) and (33,3) states of ^Hepe. At the same time, other 
characteristics of the (L — 1, A^) state of ^Hepe are close enough in comparison with the 
(L, A^) state of ^Hepe. As for Aq = 4 transitions of ^Hepe, in this case the calculated decay 
rates can be estimated only to an order of magnitude and isotope effect is uncertain. 

Thus, the Auger decay rates have been calculated for the multipolarities Aq = 3, 4. 
The method of calculation should be improved to determine with more precision the small 
components of the wave function and small decay rates in the cases of higher multipo- 
larities Ao > 4. Isotopic effect has been found and the role of the structure of the wave 
function has been studied. The Auger decay rates of metastable states Aq = 3 amount to 
tenth of the radiative transition rates and can be important in cascade processes. If the 
isotopic effect for the Aq = 3 states is the same as for Aq = 4 states, the Auger decay in 
^Hepe can compete with radiative transitions. 



7 Conclusions 

Properties of new three-body systems, the recently discovered antiprotonic helium atoms 
and analogous hadron-containing systems are intensively investigated both theoretically 
and experimentally. The investigation of intrinsic properties of antiprotonic helium atoms, 
their formation and collisions with atoms and molecules is necessary for understanding 
the antiproton fate in media. The coexistence of a particle and an antiparticle in the 
same atomic system brings on the principal possibility to study antiproton properties by 
precise measurements of antiprotonic helium atoms. 

Nonstability and an extremely large total angular momentum L ~ 30 — 40 of these 
exotic systems give rise to significant difficulties in the theoretical treatment of this system. 
However, due to very long lifetimes the metastable states of these systems are treated as 
true bound states and the variational approach is applied in the calculation. Variational 
eigenfunctions and eigenenergies are used to obtain the radiative transition rates, energy- 
level splitting due to relativistic interactions and Auger decay rates of antiprotonic helium 
atoms. The radiative transitions and Auger decay determine the cascade process after the 
formation of antiprotonic helium. Calculation of the energy-level splitting will be helpful 
in the proposed measurement of the multiplet structure [T^. One can mention that these 



rather different values are sensitive to different parts of the wave function, e. g. the region 
of small interparticle distances is important in the calculation of energy-level splitting and 
the asymptotic region is important in the calculation of Auger decay rates. 

A simple set of trial functions used in the calculations, nevertheless, provides a de- 
scription of rather different properties of antiprotonic helium atoms. Convergence of the 
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calculated values with increasing the number of trial functions has been achieved and the 
isotopic effects have been investigated by comparison of the calculated properties of ^Hepe 
and ^Hepe. In this respect, one should mention the investigation of ^Hepe having in mind 
that the reduced-mass ratio for ^Hep and '^Hep is almost the same as for '^Hep and '^Hep. 

The calculated wave functions can be used also to consider more complicated prob- 
lems; the most important of them is the determination of formation probabilities and 
initial populations. Due to substantial density and impurity effects found in experi- 
ments 0], other important applications are the interaction of antiprotonic helium 



atoms with usual helium atoms and diatomic molecules. As for the Auger decay, more 
efforts are needed to calculate small components of the wave function in the asymptotic 
region in cases of higher multipolarities Aq > 4. 
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